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Abstract. We derive a thin film model for viscoelastic liquids under strong slip which obey the stress tensor 
dynamics of corotational Jeffreys fluids. 

PACS. 83. 60. Be Viscoelasticity - 47.50.+d Non-newtonian fluid flows - 68.15.+e Liquid thin films 



Dewetting liquid polymer films on nonwetting sub- 
strates, such as silicone wafers grafted with a monolayer 
■ of brushes, play a prominent role in many nanotechno- 
logical applications. It is known that for these situations, 
polymer films on the scale of a few hundred nanometers 
typically show large slippage pp. Furthermore, for highly 
entangled polymers, the assumption of a Newtonian fluid 
will seize to be valid. To understand the interplay of vis- 
cous and viscoelastic properties of liquid polymers on hy- 
drophobically coated substrates, there is a need for refined 
theoretical methods that are able to capture and evolve 
the emerging morphologies and their longtime dynamics. 
Dimension-reduced thin film models have shown in the 
past to be extremly successful to enable quantitative pre- 
dictions that are hardly being attained simply via the un- 
derlying free-boundary problem. 

i 

In this paper we make an important step in that direc- 
tion by developing a new thin film model that combines 
large slippage with viscoelastic properties. In [2] a fam- 
ily of thin film models ranging from weak to strong slip 
regimes could be derived depending on the order of mag- 
nitude of the slip length (see also [3]). Modelling the vis- 
coelastic properties of such polymers generalized Maxwell 
and Jeffreys models have been widely used. In 0] a weak 
slip model could be combined with the linearized Jeffreys 
model to discuss effects of viscoelastic relaxation. More 
recently [H] we have shown that the strong slip limit can 
also be recovered for the linear Jeffreys model. In this 
Rapid Note we show that, for the strong slip regime, we 



are able to fully incorporate the general corotational Jef- 
freys model into our thin film model. 

We begin by presenting the underlying free boundary 
problem for incompressible, viscous flow with velocity u, 



V-u = 



du 



-Vp + V ■ r, 



(1) 



(2) 



where we assume that the traceless part of the symmetric 
stress tensor r obeys the corotational Jeffreys model |H] 



d t ( . n-i 



(3) 



Here, D/Dt denotes the Jaumann derivative which for 
arbitrary tensor fields A is given by 



DA dA 1 . A A . 

= h - (ujA - Aui) 

Dt dt 2 y ' 



(4) 



where 7 and us denote the rate of strain tensor and the 
vorticity tensor, given by 



7 = Vu+(Vu) f , u = Vu-(Vu) f , 



(5) 



respectively; d/dt is the material derivative df + u • V. We 
assume the viscosity [i as well as the relaxation parameters 
Ai, A2 to be constant material parameters. 
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As boundary conditions we assume that the substrate 
is impermeable ue z = at z — 0, and further the Navier- 
slip boundary condition 



e 2 Re*^ = -e 2 cLp + e 2 d K T xx + &T xz , (12) 



and 



u-eii = -r-en, 

1 \x " 



(6) 



with a unit- vector <l| parallel to the substrate (i.e., in x or 
y-direction) . At the liquid surface the normal component 
of the stress is balanced by the Laplace pressure (with 
surface tension coefficient 7) and the disjoining pressure 
V(z), while the tangential components of the stress tensor 
arc zero. In the following we use the reduced pressure pa = 
P+V(z). 

To simplify the algebra we restrict the calculations to 
follow to the 2D case, where we denote the velocity field 
components u = (u, w). We employ the strong-slip scaling, 
as in |2j. In this limit, the friction between the liquid and 
the substrate is too weak to maintain a non-zero xz-shear 
stress to lowest order, and lateral pressure gradients are 
balanced by the ra-component of the stress tensor. For 
the stress tensor we use the same scalings as in 0], which 
means we assume that corresponding components of the 
strain rate and stress tensor are of the same order. Hence, 
the dimensional stress tensor reads in terms of the non- 
dimensional stress tensor components 



T 



( 



X 



(7) 



where \i denotes viscosity and T denotes the time scale, 
given by T = U/L. U denotes the characteristic velocity 
scale, such as the dewetting speed and L the relative scale 
of the lateral extension of the dewetting rim. If H denotes 
the relative height of the rim then we let e = H/ 1<1. 

The dimensionless corotational Jeffreys model can then 
be written as 



1 + Al ^ 



(8) 



2//( 1 + - A 2M - e 2 {d x wf 



d z u 



dxW I T 



2m (l + A 2 |) d z w + A 2M - £ 2 (c^) 2 



(9) 



(10) 



1 + Xi^j T xz + ^(d z u - e 2 d x w){r xx - t zz ) 
= /u (l + -^2^0 (d z u + e 2 d x w) + A 2 /i (d z u - e 2 d x w) d^u . 

These equations arc coupled to the nondimcnsional gov- 
erning equations, which read in the strong slip case as 



dxU + d z w = . 



(11) 



e 2 Re*^ = -9 2 p + 4r^ + cir 2 
at 



(13) 



The corresponding scaled boundary conditions at z = 
h(x,t) are given by 



PR 



s 2 aich 2 T xx - 2d x hT xz + t z 
l + eW) 2 



{l + sHWf 12 

(14) 



and 



r xz (l-s 2 {d x h) 2 )-e 2 d x hT = Q 



(15) 



In the last equation we have introduced the first normal 
stress difference r (commonly denoted by N\ [Hj) 



Finally, the kinematic condition reads 



d t h — —V • / dz u 



(16) 



(17) 



Note that pr = p + V(h), where V(h) here denotes a 
disjoining pressure due to dispersion forces across the film. 

At the substrate, z = 0, we have as usual the imperme- 
ability condition w = and the slip condition u = br xz . 
In the strong-slip regime we have b = (3 s /e 2 , j3 s being the 
slip length parameter of order 0(e°), see j^l for details. 

For the derivation of the lubrication approximation we 
assume that the dynamic variables can be given in form 
of the asymptotic expansions 

(u, w, h,p, r y ) = (u ,w 0) h,p 0) TQ j ) (18) 
+e 2 (ui,wi,/ii,pi,T 1 J ) + 0(e 4 ) . 

To leading order in e we require from equations JSJ and 
© that 



d z u = or Tq Z = ^-d z u . 

M 



(19) 



The governing equations are, to leading order, 

dxU + d z w = , (20) 



4t^ = o, 



(21) 



d zPo = d x r* z + d z r zz , (22) 
and the leading order boundary conditions at z = ho are 

(23) 



t xz — n 



Pro ~ 2 ( ^- - 3 x h T xz ) + d xx h = , (24) 



d t ho - w + ttoci^o = . 



(25) 
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Leading order boundary conditions at the substrate are 

w = and t£ z = . (26) 

We now integrate l|21[l with respect to z and use the 
boundary conditions (|2*3|l to find 



rS z = 0. 



(27) 



With the constitutive equations l|19|) and the boundary 
conditions we are led to the plug flow condition d z uo = 
0. 

The next order equations for (JSJl , @ are then given by 



l + Ai 



2(1 + A 2 — )d z w , (29) 



where dP/dt = dt + uodc + wod z . Also, from (|22|l . we find 
with the solution (|27|l and the boundary condition l|24|) 
for the pressure at the liquid surface z — ho 

Po = r Z2 - cL x h - V(ho) ■ (30) 

The next order (0(e 2 )) it-momentum equation is 

,d°u 



Re* 



dt 



= -d xPo + d x T* x + d z Tf 



Using (|30|) and denoting uq — f(x,t) we obtain 

Re* (d t f + fd x f) = m + dci^ho + V)+ d z r{ z . (31) 

Integrating this last equation from to ho we find, using 
the slip boundary conditions to the next order, rf z = 

f/Ps, 



h Re* (d t f + fdcf) = d x / dzr - T \ z= h 4^0 (32) 
Jo 

f 

+ hoO^i^ho + V)+ r* z \ z=ho - j- . 

Ps 

The next order tangential stress boundary condition at 
z = ho is 

7f» =^^0 (33) 

Hence 



h Re* (d t f + fdxf) 



dz To 



+ hod x (8 xx ho + V)-j-. (34) 

Ps 

From equations (|28|l and (|29|) we obtain an equation for 
the difference of the diagonal terms of the stress tensor, 

(1 + Aid t + Ar/4 - Ai^/^tu" (35) 
= 4 (d x f + X 2 d t d x f + A 2 /(c*J) 2 - \2zd x fd z d x f) ■ 



We now define a film-average S of tq as 



S — -T7— / dz T . 

ih Jo 



(36) 



We denote the RHS of l|35(l by G(x,t) and observe that 
the last term is zero. Integrating IMol) with respect to z 
then yields 

Ah S + X 1 d t (4h S) + Ar/cU^oS) + X 1 (4h S)8 x f 
-\iTo~\ z =h (d t h + f&ho + hod- x f) = h G(x, t) . (37) 

Using the kinematic condition to leading order 

d t h + cLifho) = (38) 

and employing the definition of G(x, t) then finally leads 
to 

(1 + X x d t + Mf^S = (^ + \ 2 d t d x + \ 2 fd xx )f . (39) 
The lubrication model can finally be written as 

h Re* (d t f + fdxf) = d x (AhoS) + hod x (d xx ho + V)--L 

Ps 

(40) 

together with Ij38(l and l)39[l. Note that only the advective 
non-linearities, but not the corotational non-linearities ap- 
pear in (|39J) . This implies, that using the upper or lower 
advective derivative instead of the Jaumann derivative 
(leading to Oldroyd's model A or B) leads to the same 
one-dimensional thin film equation. 

We are currently investigating the effect of viscoelastic 
relaxation in strongly slipping films on the morphology of 
dewetting rims. 
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